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Abstract. - We study the phase- locking transition of two coupled low-dimensional superfluids, 
either two-dimensional superfluids at finite temperature, or one-dimensional superfluids at zero 
temperature. We find that the superfluids have a strong tendency to phase-lock. The phase-locking 
is accompanied by a sizeable increase of the transition temperature T c (in 2D systems) of the re- 
sulting double-layer superfluid to thermal Bose gas transition, compared to the Kosterlitz-Thouless 
temperature Tkt of the uncoupled 2D systems, which suggests a plausible way of observing the 
Kibble-Zurek mechanism in two-dimensional cold atom systems by rapidly varying the tunneling 
rate between the superfluids. If there is also interaction between atoms in different layers present 
we find additional phases, while no sliding phase, characterized by order or quasi long range order 
(QLRO) either in the symmetric or the antisymmetric sector of the system. 



Introduction. — Since the realization of the Mott 
insulator transition [1], remarkable progress has been 
made in controlling and manipulating ensembles of ultra- 
cold atoms [2], which was followed by a number of ex- 
periments to create and study more and more sophis- 
ticated many-body effects, such as fermionic superfluids 
. [3], one-dimensional strongly correlated Fermi and Bose 
systems [4], or noise correlations in interacting atomic 
systems [5, 6]. An intriguing new direction in study- 
ing low-dimensional strongly correlated systems was taken 
with the realization and observation [7] of the Kosterlitz- 
Thouless (KT) transition [8]. In this experiment the 
interference amplitude between two independent two- 
dimensional (2D) Bose systems was studied as a function 
of temperature. This analysis revealed the jump in the su- 
perfluid stiffness (see also Ref. [9]) and the emergence of 
unpaired isolated vortices as they crossed the phase tran- 
sition. 

In another interesting recent experiment Sadler et. al. 
observed spontaneous generation of topological defects in 
the spinor condensate after a sudden quench (i.e. a rapid, 
non-adiabatic ramp) through a quantum phase transi- 
tion [10]. A similar experiment in a double-layer system 
was reported in Ref. [11]. The topological defects are gen- 
erated [12] at a density which is related to the rate at 
which the transition is crossed [13]. Later it was argued 
that the dependence of the number of such defects on the 



swipe rate across a quantum critical point can be used as 
a probe of the critical exponents characterizing the phase 
transition [14]. This Kibble-Zurek (KZ) mechanism was 
originally considered as an early universe scenario creating 
cosmic strings, which would serve as an ingredient for the 
formation of galaxies [15]. Cold atom systems appear to 
be a very suitable laboratory for performing such "cos- 
mological experiments" , since these systems are highly 
tunable and well isolated from the environment. So far 
the experiments and the theoretical proposals addressed 
the KZ scenario across a quantum phase transition. The 
main reason is that it is generally hard to cool such sys- 
tems sufficiently fast to observe non-equilibrium effects. 
In this work we provide an example of a particular system 
where this difficulty can be easily overcome by quench- 
ing the transition temperature T c instead of T. Thus the 
relevant ratio T/T c can be tuned with an arbitrary rate 
and the KZ mechanism can be observed. Specifically, we 
examine a system of two superfluids (SF): As we show be- 
low, by turning on tunneling between the two systems the 
transition temperature increases rapidly, and the system 
attempts to create long-range order (LRO). However, in 
this process, defects in the SF phase are created, which 
develop into long-lived vortex-anti-vortex pairs or in fi- 
nite system unbalanced population between vortices and 
anti- vortices. We note that because the systems are iso- 
lated and there is no external heat bath, the temperature 
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Fig. 1: Phase diagrams of two 2D SFs, coupled through a 
term of the form S12, Eq. (5), in terms of Jint/J and T/Tkt- 
For low temperatures we find antisymmetric quasi-order (AQ) 
and/or symmetric quasi-order (SQ), which either simultane- 
ously undergo a KT transition due to single vortices (AQ/SQ 
to thermal Bose gas (TBG) phase), or individually due to cor- 
related vortex pairs: symmetric (anti-symmetric) vortex pairs 
drive the AQ/SQ to AQ (SQ) transition. 

itself also changes due to the quench. However, the long- 
wavelcngth fluctuations relevant for the KT transition are 
only a small subset of all degrees of freedom, majority of 
which are only weakly affected by small inter-layer tunnel- 
ing. So we believe that the change of the T c is the main 
effect of the quench. 

In this work we consider two SFs coupled via tunnel- 
ing and/or interactions. In the experiments the hopping 
or tunneling rate between two systems can be tuned to 
a high precision [4,7,16,17]. Interactions between the 
atoms in different systems can either be realized in en- 
sembles of polar molecules or by using mixtures of two 
hyperfine states, where the tunneling rate is controlled by 
an infrared light source [18], which induces spin- flipping 
between the hyperfine states. In this case the atoms in dif- 
ferent states naturally interact with each other since they 
are not physically separated in space. We also discuss the 
analogous one-dimensional (ID), zero temperature lattice 
systems. The main results of our analysis are the phase 
diagrams of coupled SFs in Fig. 2 (2D) and Fig. 4 (ID), 
the behavior of T c and the energy gap shown in Fig. 3, 
as well as the proposal of realizing the KZ mechanism by 
switching on the tunneling between two SFs. 

2D superfluids. — In this section we consider two 
2D SFs, each characterized by a KT temperature Tkt- 
We write the bosonic operators b 1 / 2 in the two lay- 
ers in a phase-density representation [8, 19], b\/ 2 ~ 
y/pi/2 ex P(*0i/2)) where pxji are the density operators of 
the two systems, and (f>i/ 2 the phases. The low-momentum 
fluctuations of the phase fields are described by Gaussian 
contributions to the Hamiltonian Ho- Because of the for- 
mal analogy between the quantum ID and thermal 2D 
systems [23] we adopt the quantum terminology through- 
out the paper and refer to the ratio of the Hamiltonian 
and the temperature as the action. Then 

* = ^ = ^/*r[(Vfc) 9 + (Vfc) a ]. (1) 



The energy scale J here is related to Tkt by J = 2Tkt/k- 
Besides these long-wavelength fluctuations, the system 
also contains additional degrees of freedom, vortex-anti- 
vortex pairs [8]. The corresponding term in the action is 
expressed through the dual fields 9 1,2 [19]: 

Si = (2^ [cos(2 ^ )+cosM ' (2) 

where a is a short-distance cut-off of the size of the vor- 
tex core, and A\ is proportional to the single-vortex fu- 
gacity: A\ ~ Jexp(— J/T), where we assume both SFs 
to have the same effective parameters J and A\. Op- 
erators of the type exp(2i9) create kinks in the field <j>: 
&cg,{-2i9{x))(j){x')ex^{2i9{x)) ~ 4>{x')+2-kQ{x-x'), &{x) 
being the step function, which corresponds to the effect of 
vortices in the original 2D problem (Ref. [23], p. 92). 

In addition the two systems are coupled by a hopping 
term ~ t±b\b 2 + h.c, which results in the following con- 
tribution to the action: 

S ^ = W(2^ C ° S( ^- 02) ' (3) 

where the bare value of J± corresponds approximately 
to t±po- In principle, the hopping term is modified by 
the vortex contributions, however, these corrections are 
always irrelevant under rcnormalization group (RG). 

For most of the discussion in this paper we use the sym- 
metric and anti-symmetric combinations of <p 1 / 2 and 9 1 / 2 - 

<P s/a = (<j> 1 ± <h)/y/2, 9 s/a = (9 1 ±9 2 )/V2. (4) 

Written in these fields, the term So in Eq. (1) is again 
a sum of Gaussian models, now in the fields (j) s and <f> a , 
with the same energy scale J. However, we will consider a 
broader class of actions, in which the energy scales of the 
symmetric and anti-symmetric sector differ. We include 
the following term in the action: 

S12 = ^ J d 2 rV0!V0 2 (5) 

With this, the quadratic part of the action is given by: 

So+S 12 = A| d2r (V^) 2 + ^|rf 2 r(V0 a ) 2 ,(6) 

where J 8 and J a are given by J s / a = J ± J ln t- 

We now motivate the existence of such a term S12 in 
ultracold atom systems, by considering two BECs cou- 
pled by a short-range density-density interaction. Start- 
ing from a Hamiltonian of the form H = X)fe[ ek ^k^ k + 
(g/2V)p^pk\, where 6k is the boson operator, ek the free 
dispersion ek = k 2 /2m, g is the interaction strength of the 
contact interaction, V the volume, and pk is the density 
operator of momentum k, given by p^ = ^ p fe^frp+k, we 
assume that the zero momentum mode is macroscopically 
occupied, and formally replace the operator bo by a num- 
ber, bo — » ^/No, where N Q is the number of condensed 
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Fig. 2: Phase diagram, temperature (in units of Tkt) ver- 
sus interaction (in units of J). We assume J±/J ~ 1CP 3 and 
A x I J ~ f(T 3 . DLSF: double layer superfluid; TBG: ther- 
mal Bose gas; ASF: anti-symmetric superfluid; SQ: symmetric 
quasi-order. The order of the transition lines are either first 
(I), second order (II), or KT (thin lines). 



atoms which is comparable to the total atom number N, 
i.e. Nq < TV. Next we keep all terms that are quadratic 
in fek (with k ^ 0), and perform a Bogoliubov transfor- 
mation, given by: &k = «k/3k + Wk/?L k , to diagonalize the 
Hamiltonian. The cigenmodes /?k have a dispersion rela- 
tion = y/ £k(ek + 2gn), with n being the density N/V. 
The low-k limit is given by ~ w 2 |k| 2 , with v = \J gn/m, 
which corresponds to the contribution in Eq. (1) of the 
action. Next, we consider the sum of two copies of the 
previous Hamiltonian with boson operators b 1 / 2 - In addi- 
tion we consider an interaction H\2 = 9i2/VJ2kPi k/°2,k, 
where the density operators pi/2,k are given by Pi/2,k = 
S P &1/2 p^i/2.p+k- Following the same procedure as be- 
fore, we find two eigenmode branches, corresponding to 
in-phase and out-of-phase superpositions of the modes of 

kl 2 , 



each condensate, with the dispersions w 2 / a k 



J s/a\ 



with the velocities v s / a = y/(g ± gi2)n/m. Therefore, for 
this example, the energy scale Jj„ t is related to gi2n/m, 
which would be of similar order as J for a system interact- 
ing via contact interaction, for small temperatures. This 
discussion only applies to the weakly interacting limit of a 
true condensate. However, it demonstrate that a density- 
density contact interaction term can lead to a substantial 
energy splitting of the in-phase and out-of phase modes. 

Finally, in addition to single vortices in each SF, we have 
to consider the possibility of correlated vortex pairs, i.e. 
one vortex in each layer at the same location of either the 
same or of opposite vorticity. We will refer to these vor- 
tex configurations as symmetric or anti- symmetric vortex 
pairs, respectively. These excitations appear as the follow- 
ing terms in the action: 



S s a — 



2A S 



d 2 r 
(27ra) : 



cos(2\/2# s , Q ). (7) 



the action, which derives from the possibility that the vor- 
tices in different layers interact with each other, through 
the terms (3) and (5). The effect of these terms is the fol- 
lowing: At low temperatures the energy between two sin- 
gle vortices of opposite vorticity due to tunneling grows 
as the square of the distance D between them, i.e. as 
J±(D/a) 2 . As a result, the tunneling term attempts to 
confine vortices of opposite vorticity, leading to phase- 
locking between the layers, which we describe further later 
on. The interaction S12 changes the energy of correlated 
vortex pairs as follows: The energy of a single vortex is 
given by 27rJlogL/a, where L is the system size, whereas 
a symmetric/anti-symmetric vortex pair has an energy of 
47r( J ± Ji n t) logL/a. Therefore, symmetric vortex pairs 
are the lowest energy vortex excitations for Jj nt < — J/2, 
whereas for J int > J/2 anti-symmetric vortex pairs are the 
lowest energy excitations. As we will see below, in these 
regimes correlated vortex pairs drive transitions to phases, 
in which one sector is (quasi-)SF whereas the other is dis- 
ordered. We will also see that these terms are generated 
under the RG flow, even if not present at the onset. 

We note that a similar system has been studied in [28]. 
Here, we consider a larger class of systems by including the 
interaction term (5), which in turn requires us to include 
the correlated vortex excitations (7). These terms give 
rise to additional phases as we will see in the following. 

Next we analyze our system within the RG approach. 
This RG flow is perturbative in the vortex fugacities A\, 
A s , and A ai and the tunneling energy Jj_, and therefore 
applies to the weak-coupling limit (in particular J± — ► 
+0). At second order the flow equations are given by [21]: 



M~( 2 2nJ a ) J 



(8) 



These correlated vortex terms, which describe new degrees 
of freedom, can be the most relevant non-linear terms in 



— - (2 - 2n- j A s + , (9) 

dA a ( J a \ Aj(J s - Jq) 

— = [2 2*- j A a + a 3 ^ 2 , (10) 

^ = -a 2 (2^J 2 + ^(J s + J a )J s )2J s . (13) 

The coefficients a 2 /3 are non-universal parameters that 
appear in the RG procedure [22], and which do not af- 
fect the results qualitatively. For consistency, we have to 
expand the right-hand site of the above equations up to 
second order, around the resulting Gaussian fixed point: 
J si a — J ± Jint + is/a- We emphasize again that Jj„ t near 
the fixed point can be generated by RG and be nonzero 
even if it is not present at the onset. 

Before we consider the full RG flow, we consider the sim- 
pler case of no tunneling, i.e. we solve the RG equations 
while setting J± = 0. In Fig. 1 we show the phase diagram 
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of two 2D SFs coupled by S12, Eq. (5). Such a system 
would be realized by a 2D mixture of bosonic atoms in 
two different hyperfme states, interacting via some short- 
range potential. The order parameters we consider are 
O s (x) = 6i(x)6 2 (a ; ) and O a {x) = b\(x)b2(x). To obtain 
the phase diagram we consider the correlation functions 
of each of these order parameters, which can either scale 
algebraically or exponentially. In Fig. 1 we refer to alge- 
braic scaling of O s (x) as symmetric quasi-order (SQ), and 
of O a (x) as anti-symmetric quasi-order (AQ). In each of 
the sectors a KT transition marks the transition from the 
algebraic to the exponential regime, which occur either si- 
multaneously and are driven by single- vortex excitations, 
or at different temperatures and are driven by correlated 
vortex pairs. As a result we find four regimes: At temper- 
atures above Tkt, both sectors are disordered, giving rise 
to a thermal Bose gas (TBG) phase. For temperatures be- 
low Tkt, and for a wide range of Ji n t, we find that both 
sectors are quasi SF (AQ/SQ), which is the only phase in 
which the correlation function of the single boson opera- 
tors show algebraic scaling. We also find regimes in which 
only one sector shows algebraic scaling, whereas the other 
is disordered (AQ and SQ). From the perspective of vor- 
tices, the TBG phase is a gas of free single vortices in each 
layer, whereas the AQ (SQ) phase is a gas of symmetric 
(anti-symmetric) vortex pairs. 

We now consider the full RG system, including J±. 
We numerically integrate the RG equations, and find the 
phase diagram shown in Fig. 2 in terms of the tempera- 
ture T and the interaction J; nt . We again find four differ- 
ent phases that are different combinations of LRO, QLRO, 
and disorder in the symmetric and anti-symmetric sector. 
At high temperatures we find that both sectors are disor- 
dered in a TBG phase, as before. For lower temperatures, 
and for a wide range of J; n t , the system is in a double-layer 
SF phase (DLSF): The symmetric sector shows algebraic 
scaling, whereas the exponent of the anti-symmetric sector 
is renormalized to zero, i.e. we find two SFs that are phase- 
locked due to J±. Note that the transition temperature 
T c between DLSF and TBG has been noticeably increased 
relative to the decoupled value Tkt, as we will discuss fur- 
ther later on. We also find two additional phases, which 
are partially (quasi-)SF and partially disordered. One of 
them is the SQ phase, as before, whereas the other one 
(ASF), now shows true LRO in the anti-symmetric sector 
due to J±, whereas the symmetric sector remains disor- 
dered. We note that the generic double-layer action that 
we discuss in this paper does not show a sliding phase [20] , 
for any non-zero Jj_. Either Si or S a , which is generated 
by RG, drives the anti-symmetric sector to a disordered 
state, or S± creates true LRO in the field (j> a . 

We also use the RG flow to find the order of the 
phase transitions in the weak-coupling limit that the anti- 
symmetric sector undergoes, by determining the energy 
gap using a 'poor-man's scaling' argument: when the cou- 
pling amplitude J±(l*) is of order unity the corresponding 
gap is given by the expression A ~ J±exp(—l*). From 
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Fig. 3: (a) Critical temperature T c of the DLSF- TBG transi- 
tion (in units of Tkt) for different values of A\/J: 10~ 3 , 0.1, 0.4 
(I— III) , and for Jj n t = 0. (b) Energy gap in the anti-symmetric 
sector (in units of J±) as a function of J±/J and temperature 
(in units of Tkt)- We have set Ai/J — 0.1 and Jint = 0. 

the behavior of A at the phase transition we can read off 
whether it is of first or second order, as indicated in Fig. 2. 

Given the nature of an effective theory, only approx- 
imate statements can be made about how the different 
regimes of the phase diagram relate to the microscopic in- 
teractions. To create the ASF or AQ phase an attraction 
between the two atom species is needed that is of order J, 
whereas to create the SQ phase, a repulsion of that order 
would be needed. To detect the different phases, one could 
use the interference method used in [7] to distinguish the 
phase-locked phases (DLSF and ASF), which would show 
a well-defined interference pattern, from the uncorrelated 
phases. Another approach would be time-of-flight images: 
The DLSF phase would display a quasi-condensate signa- 
ture, whereas the other phases would appear disordered. 
However, at the transition from ASF or SQ to TBG, the 
width of the distribution would abruptly increase. 

Kibble-Zurek mechanism. — In this section we dis- 
cuss how the phase-locking transition found in the previ- 
ous section could be used to realize the KZ mechanism. 
The defining property of this mechanism is the generation 
of topological defects by ramping across a phase transition, 
coming from the disordered phase. The disordered phase 
that we propose to use is the TBG phase of the decoupled 
2D systems, that is, we consider the experimental setup 
reported in ref. [7] for a temperature T above the KT tem- 
perature Tkt- The ordered phase we consider is the DLSF 
phase, i.e. the phase-locked phase of two coupled SFs. The 
ramping is achieved by turning on the tunneling between 
the two layers, which can be done by lowering the poten- 
tial barrier between them. For this procedure the critical 
temperature T c of the DLSF-TBG transition needs to be 
above the KT temperature of the uncoupled systems. We 
now show that the RG flow indeed predicts such a scenario. 
In the experiments in Ref. [7] , the atoms in different layers 
do not interact with each other. Therefore, it can be ex- 
pected that Jint is small, of order Jj_, which motivates us 
to discuss the case Ji nt = here. We note however, that 
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the desired scenario of an increased critical temperature, 
is found for a wide range of Jmt, as can be seen in Fig. 2. 
In Fig. 3 (a) we show how the critical temperature of the 
DLSF-TBG transition behaves, predicted by the RG flow, 
for different values of A\. The critical temperature shows 
a sizeable increase, due to the phase-locking transition. 
Due to the perturbative nature of the RG scheme, the RG 
flow underestimates the effects of the term S±, and pre- 
dicts a finite jump of the critical temperature when Jj_ is 
turned on. However, to lock the SFs together in the regime 
slightly above Tkt, J± needs to be at least of the order 
of the vortex core energy, giving rise to a finite slope of 
T c instead of a jump. The energy gap of this transition is 
shown in Fig. 3 from which we can see that the transition 
is of first order, in contrast to the second order transi- 
tion described in [12, 13], which is advantageous because 
the onset of order is instantaneous rather than continu- 
ous. We note that the phase diagram was obtained using 
the assumption that the bare parameters of the model, in 
particular J, do not depend on temperature. This is true 
only if temperature is close to Tkt- Here we find that the 
ratio T c /Tkt can be relatively large. In fact T c /Tkt will 
be always smaller than that shown in Fig 3 (a), however, 
qualitatively the behavior of T c /Tkt as a function of J± 
should remain intact. We point out that our results can 
be generalized to a system of N > 2 coupled SFs. One 
finds that the SFs still show a strong tendency to phase- 
lock together. As a result the critical temperature should 
approximately satisfy the equality ttJ(T c )N — 2T C . Thus 
as N increases T c approaches the mean-field critical tem- 
perature at which the stiffness J vanishes and we recover 
the usual 3D result. 

In finite size systems there is another constraint on the 
minimum value of J± : We consider the free energy of a sin- 
gle vortex in the anti-symmetric field: <j) a ~ arctan(cc/y). 
For the decoupled system we get for the free energy [22]: 
F <~ 2(nJ — 2T)\ogL/a 1 where L is the system size. 
The coupling term gives a free energy contribution F± ~ 
J±(L/a) 2 . In the thermodynamic limit, L — ► oo, this 
term diverges faster than the others, which is consistent 
with our finding of LRO in the antisymmetric sector. For 
a finite system, comparing these terms gives the estimate 
Jj_ <~ J log(i / a) I (L I 'a) 2 , that is required for this order 
to develop. With a system size L/a > 10 2 , that would 
require J± > 10 _3 J, which, for the setup in [7], would be 
around 10 2 s _1 . 

As an estimate of the number of domains that would 
be created, we follow the argument in [12]: The coher- 
ence scale of the DLSF phase is given by (A/J) 1 / 2 ^, 
which is the scale of a Klein-Gordon model with a ki- 
netic energy scale J and a 'mass-term' with a prefactor 
A/a 2 . The domain size is then given by (A/ J) a 2 , and 
the number of domains by ~ (J/A)L 2 /a 2 . As we show 
in Fig. 3 b) for J ± /J w KT 2 , we find A/J± ~ KT 1 , 
and therefore J/ A ~ 10 3 . With L/a > 10 2 , we would get 
Ndom ~ 10 1 — 10 2 , which would generate a similar number 
of vortices. We estimate the vortex-antivortex imbalance 
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Fig. 4: Phase diagram of two coupled lattice ID bosonic SFs at 
unit filling: Luttinger parameter K versus interaction (in units 
of SF velocity). MI: Mott insulator; SQ: symmetric quasi- 
order; PS: phase separation; CSF: coupled superfluids; CL: 
collapsed phase; ASF: anti-symmetric SF. 

by considering the number of domains around the periph- 
ery of the system, which scales as L/£. If we imagine that 
the phase behaves like a random walk, the total phase mis- 
match, corresponding to the vortex-antivortex imbalance, 
will scale as \f~L~J\ ~ N]^, which, for L/a > 10 2 , is of 
the order 10° - 10 1 . 

In summary, we propose the following procedure: i) Pre- 
pare two uncoupled SFs at a temperature T slightly above 
Tkt- h) Switch on the tunneling between the two layers, 
which creates a DLSF phase with a critical temperature 
T c higher than T. As a result, one should find a number 
of long-lived vortex-antivortex pairs in the anti-symmetric 
phase field <f> a , which would be visible in an interference 
measurement, at a temperature where there would be none 
in thermal equilibrium. 

ID superfluids. — We now study the phases of two 
coupled SFs in ID (see [25-27]). We consider both a lattice 
system with unit filling, and a system with incommensu- 
rate filling. These types of systems are formally closely 
related to the 2D system that we discussed before [23]. 
The imaginary time direction in a quantum ID system 
plays the role of the second dimension in the 2D case. 

The action of two incommensurate SFs, coupled by a 
hopping term, is given by S = So.i + So, 2 + S±. In 
the terms 50,1/2 we formally replace the ratio irJ/T by 
the Luttinger parameter K, and similarly we introduce 
g± = J±/T in the hopping term. The parameter K and 
a velocity v, which is contained in r = (vt,x), charac- 
terize these ID SFs, where if is a measure of the in- 
teraction: weakly interacting bosons have K » 1, while 
bosons interacting via a strongly repulsive potential, have 
K — ► 1 (Ref. [28]). The system again separates into sym- 
metric and anti-symmetric fields, Eq. (4), where the anti- 
symmetric sector is given by a sine-Gordon model. By 
going to this basis, one can determine that the scaling di- 
mension of the operator cos(\/20 a ) is given by 1/(2K). 
Therefore for any value of K, this term is highly relevant, 
and the anti-symmetric sector develops LRO, K a — > oo, 
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giving rise to a phase of coupled SFs (CSF). 

For a lattice system of SFs with unit filling, this phase- 
locking tendency competes with the Mott insulator tran- 
sition, and other, partial, localization transitions. We 
include the terms S± and S a a in the action, where we 
replace A\jT by g\, and A s ^ a /T by g s ^ a , and finally we 
add S[ 2 = (Ui2/n 2 ) J d 2 r\79iV9 2 which would derive from 
some inter-chain density-density interaction. From these 
expressions it is clear that the ID commensurate system is 
formally identical to the 2D system at finite temperature. 
We again integrate the RG equations numerically, and we 
find that if the parameters g± and gi are of comparable 
magnitude then the hopping term is the most relevant 
one, driving the system to a CSF phase. However, for 
<7i ^ ff-Li an d f° r small K, the terms Si,S s , or S a can 
dominate, and localize either the symmetric or the anti- 
symmetric sector, or both. The resulting phase diagram 
is very similar to that of the 2D system with the localized 
phases corresponding to the thermal ones (see Fig. 4). 

Conclusion. — In summary, we studied the phase- 
locking transition of ID and 2D SFs, within an RG ap- 
proach. For 2D SFs, this transition is accompanied by an 
increase of the transition temperature. We suggest that 
this effect can be used to probe the KZ mechanism in 
cold atom systems by rapidly changing the ratio T c /T. 
For both 2D and commensurate ID systems, interaction 
terms give rise to additional phases, in which either the 
symmetric or the anti-symmetric sector is disordered or 
localized, respectively. 
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